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I study the dynamics of a superfiuid vortex in a random potential, as in the inner crust of 
a neutron star. Below a critical flow velocity of the ambient superfiuid, a vortex is effectively 
immobilized by lattice forces even in the limit of zero dissipation. Low-velocity, translatory motion 
is not dynamically possible, a result with important implications for understanding neutron star 
precession and the dynamical properties of superfiuid nuclear matter. 
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A neutron star (NS) is expected to comprise over a so- 
lar mass of distinct quantum liquids. In the inner crust 
(between the drip density and approximately half nuclear 
density), a JS* neutron superfiuid (SF) threaded by an ar- 
ray of quantized vortices coexists with the ionic lattice. 
Forces exerted by the vortex array on the lattice could 
have important observable effects on NS spin and ther- 
mal evolution (e.g., [TJ |5J [3J |U [5]), especially if vortices 
pin to the lattice as suggested long ago [B]. Evidence 
for long-period precession in some NSs (e.g., [7]), an in- 
terpretation supported by quantitative modeling [8], is 
difficult to explain if vortices pin in the inner crust [2] 
and presents challenges to the standard picture of the 
outer core [3]. The critical question of whether or not 
pinning occurs in the inner crust has been studied (e.g., 
[Jill HI]) but not yet satisfactorily answered. Here I show 
that a vortex is trapped by the inner-crust lattice under 
very general conditions if the ambient SF velocity is be- 
low a critical value. The problem of vortex dynamics and 
pinning in a lattice potential is also of considerable in- 
terest in laboratory Bose-Einstein condensates (see, e.g., 

The vortex drag description-Most hydrodynamic 
studies of the coupling problem begin by including a term 
in the SF acceleration equation for the "mutual friction" 
force in a homogeneous medium. For a SF of mass den- 
sity p s flowing at velocity v s , coupled to a medium mov- 
ing at velocity v m , and neglecting the collective tension 
of the vortex lattice, the mutual friction force density in 
the non-rotating frame is |13l 114] 

F m f = ft p s UJX.(v s -V m )+Pp s Vx[u} X (v s - V m )] , (1) 

where oj = V x v s is the SF vorticity, u is a unit vector 
in the direction of u>, and (3 and /3' are coefficients to 
be determined by a microscopic calculation for the appli- 
cation of interest. This drag term, or extensions of it to 
bulk neutron matter in beta equilibrium, has been widely 
used in studies of NS hydrodynamics and precession; the 
"medium" could be taken to be, for example, the lattice 
of the NS inner crust, magnetic flux tubes of the outer 
core with which the vortices interact, and the charged 
fluid of the outer core. The mutual friction force of eq. 



[I] is directly related to the drag force per unit length on 
a vortex [JJ 

fd = ~V V - V' u x u > ( 2 ) 

where v is the vortex velocity with respect to the 
medium, and the coefficients ry and rf are related to (3 and 
(3' . The second term in eq. [2] is non-dissipative and is 
usually assumed to be zero, as I also assume; in this case 
(3 = T] r /(l + rft) and /?' = rj r (3, where r\ r = r]/p s K is the 
reduced drag coefficient and k is the vorticity quantum 
h/2m n (m n is the neutron mass for a neutron SF). Dissi- 
pative processes that act on a translating vortex include 
electron scattering jTSJ [TB] , and the excitation of vortex 
waves (Kelvin modes) through the vortex-nucleus inter- 
action [TTJ [TB] . All of these calculations give r\ r < < 1 . 

Neglecting vortex bending and local non-dissipative 
forces, the motion of a vortex segment follows from equat- 
ing the sum of the Magnus and drag forces to zero: 

p s K x (v - v s ) - r\v = 0, (3) 

where k is aligned with the vortex (the z axis) and v s is 
now the velocity of the ambient SF in the rest frame of 
the medium, taken to be along the —x axis. In this de- 
scription (henceforth, the "drag description"), the vortex 
segment moves at velocity v ~ v s for low drag (r/ r << 1) 
and v ~ for high drag (rj r >> 1). The large drag 
limit has been taken as corresponding to the pinning 
of vortices to defects, such as nuclei in the inner crust 
or magnetic flux tubes in the outer core [TJ [19j |2Ql EJJ. 
The drag description embodied in eq. kj], however, in- 
correctly describes the motion of a vortex through these 
fundamentally inhomogeneous environments. A correct 
description of vortex motion in the inner crust requires 
the inclusion of two additional forces: 1) the local, non- 
dissipative component of the force exerted on the vor- 
tex by the lattice, and, 2) the elastic force of the vortex. 
These forces qualitatively change the vortex motion. Ap- 
propriate to the inner crust (but not the core) , I focus on 
a single-component SF. 

Vortex motion through a random lattice-Consider the 
motion of a vortex that is dissipatively coupled to a back- 
ground lattice. The displacement vector of the vortex 
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with respect to the z-axis is u(z, t) = u x (z, t)x+u y (z, t)y. 
The equations of motion for a vortex moving in the ab- 
sence of external forces can be found in Sonin |22j : a 
hydrodynamic description suffices for excitation wave- 
lengths significantly larger than the SF coherence length 
£, about 10 fm in the inner crust. The force per unit 
length exerted on the vortex by the lattice has a non- 
dissipative contribution fo and a dissipative contribution 
taken here to be the drag force of eq. [2] with 77' = 0, 
assumed to hold locally, and approximated as linear in 
the local vortex velocity (but easily generalized). The 
equations of motion become 

m d 2 u ( du \ „ du 

T^+Ps^[^-Vs)+fo-V^=0, (4) 

where T v = (p s K 2 /47r) ln^fc)^ 1 is the vortex self-energy 
(tension) for an excitation of wave number k. The first 
term represents the restoring force of tension as the vor- 
tex is bent. The second term is the Magnus force per unit 
length exerted on a vortex that is moving with respect to 
the ambient SF. In the absence of external forces, the vor- 
tex would remain straight and motionless with respect to 
the solid. An essential feature of vortex dynamics is that 
the local vortex velocity is determined entirely by exter- 
nal forces and by the shape of the vortex, and not by 
inertial forces (for temperatures small compared to the 
pairing gap). In the absence of external forces, the solu- 
tions to eqs. [4] are circularly-polarized, diffusive waves 
(Kelvin waves) with frequencies uj^ — ±T v k 2 / p s n |23j . 
Characteristic bending scales of the vortex in the prob- 
lem are ~ 500 fm (see below), and fixing In(^k)^ 1 = 3 is 
an adequate approximation. I ignore quantum effects on 
the vortex motion and thermal excitations. 

Nuclei in the solid with which the superfluid coexists 
exert local forces on the vortex which vary over a length 
scale of order the average nuclear spacing a, typically 
30-50 fm in the inner crust. Calculations of the vortex- 
nucleus interaction energy E vn give values of up to ~ 1 
MeV, attractive in some regions and repulsive in others 
[2~il |2"5] ; the corresponding characteristic force on a vor- 
tex segment of length a is F m = a~ 1 E vn . The solid is 
most likely amorphous 26J, in which case the vortex in- 
teracts with an effectively random potential Even if the 
lattice does possess long-range crystalline order, vortices 
will be aligned with the SF rotation axis on average which 
will not, in general, coincide with any symmetry axis of 
the solid; a random lattice approximation should be ad- 
equate also in this situation. In the denser regions of the 
NS inner crust, T v is of order an MeV fm -1 . Defining a 
dimensionless strength parameter s = F m /T v (> for an 
attractive potential) , estimates of F m from Refs. [2~il |2"5"] 
give 10 -2 < |s| < 0.1. A vortex is thus a stiff object; 
the vortex-nucleus interaction can bend a vortex by an 
amount a only over a length scale of order > 10a [27] ; a 
vortex segment of length a can be regarded as straight to 
a good approximation and the lattice force on the vortex 
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FIG. 1: Left panel: vortex position, averaged over its length, 
in the presence of an external SF flow field of velocity v 3 = 
10 -3 and no dissipation over 5 x 10 3 time units. Without the 
random potential, the vortex would translate to the left with 
the SF by five units. The non-dissipative lattice force effec- 
tively pins the vortex, confining it to a meandering trajectory 
over a length scale of order a. Right panel: Length-averaged 
vortex speed in response to local SF flow at v s . Velocities 
have been scaled by s~ 3 ^ 2 . For relatively high v s , the vortex 
translates at a speed ~ v s , but is trapped for v a < v c ~ s 3 ^ 2 , 
independent of drag. 

can be taken as acting in the x — y plane at any given z 
along the vortex. For the force exerted by a single nu- 
cleus at the origin on a vortex segment of length a, I take 
a parameterized central force in the x — y plane: 

f vn (u) = -1.7F m -e-" 2 ^ (5) 

The force has a maximum magnitude of F m at u = r p , 
where r p is the effective range of the potential, henceforth 
taken to be r p — a/2. Dividing the vortex into segments 
of length a, let the total non-dissipative force on the i-th 
segment be 

fo(z%) = f™( u i ~ r i)' ( 6 ) 

3 

where the nuclei are randomly placed at locations Tj in 
planes separated by a and parallel to the x — y plane. For 
simplicity, each summation is only over nuclei in the i-th 
plane. In the following, lengths will be expressed in units 
of a, time in units of tjy = p s na 2 JT V (the characteristic 
diffusion time of a Kelvin wave of wavenumber k = aT 1 
along the segment) and velocities in units of a/tjj. As 
described below, the sign of s is unimportant for r p ~ 
a; I choose s > 0, corresponding to attractive nuclear 
potentials, for illustration. I focus on the regime 0.1 > 
rj r > which widely brackets the range found by all drag 
calculations, and take s — 0.1. 

Comparison of terms in eq. [4] shows that the dynam- 
ics are essentially determined by the non-dissipative lat- 
tice force and s in the regime v s < v c . To determine the 
character of the dynamical regimes of vortex motion and 
the scaling relations that define them, eqs. [4]- [6] were 
solved numerically. The vortex was divided into 50 zones 
of length a, with periodic boundary conditions applied at 
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FIG. 2: Motion of an example vortex segment through one plane of the lattice for rj r = 0.1. Left panel: For v a > v c ~ s 3//2 , 
the vortex translates with the superfluid, with significant changes in direction as it interacts with nuclei (denoted by circles). 
For v s < v c , the segment damps to a pinned configuration that brings it closer to nuclei. Middle panel: Comparison of the 
damped motion for v 3 — and v s < v c , showing that finite v s has little effect on the motion in this regime. Right panel: Detail 
of the box region of the middle figure. Circles on the lines indicate points of equal time. The result of finite v s is to displace 
the pinned vortex slightly in the y direction. 



the ends, in a random lattice consisting of 10 3 nuclei per 
zone. Fig. [I] (left panel) shows that the vortex is trapped 
by the random potential for finite v s even for zero drag. 
The chief result of the numerical analysis is shown in the 
right panel of Fig. 1. Trapping occurs for v s below a crit- 
ical velocity v c ~ s 3 / 2 , independent of drag in the regime 
7] r << 1. For v s > v c , the vortex moves at v ~ v s ; in 
this regime, the drag description of eq. [3] is adequate. 
Below v c , however, the vortex velocity is effectively zero. 
The drag description fails to show this transition to the 
pinned state because it excludes both the non-dissipative 
lattice force and the finite vortex tension. The value of 
the critical velocity v c ~ s 3 / 2 follows from consideration 
of a vortex in static equilibrium, as discussed elsewhere 
[TO]. 127] , Only if the vortex had infinite tension, one of 
the assumptions implicit in the drag description and eq. 
|3j, would translatory states exist for any v s . 

Dissipative processes will damp the vortex to a station- 
ary, pinned configuration. Fig. 2 shows an example seg- 
ment of the long vortex moving under drag. For v s < v c , 
the trajectory is essentially the same as for v s = since 
the non-dissipative force dominates the global Magnus 
force in this velocity regime. Finite v s (below v c ), pro- 
duces only a small displacement of the vortex segment's 
final, damped position. Inspection of eq. [4] suggests 
that a dragged vortex will damp to a pinned position 
over a characteristic time ~ (sr) r )~ , and numerical ex- 
periments confirm this. The damping time is ~ 100 for 
the examples shown in Fig. 2. Fig. 3 shows the initial 
motion and damping of the vortex in the x — z plane. 
The initially-straight vortex forms kinks immediately in 
the presence of the random potential, and by t = 10 the 
vortex has been excited to amplitudes of order a, the 
length scale of the random potential. By t = 10 3 , the 
vortex has damped to a stationary pinned configuration 
with bends over a characteristic length scale of ~ 10a. 
Because the vortex has large tension, it cannot bend to 
intersect every nucleus, but assumes a shape that strikes 



the best compromise between the energy gain of being 
close to a nucleus and the cost of bending the vortex. 

Further simulations show that these conclusions are 
unaffected if the attractive nuclear potentials are re- 
placed by repulsive ones (though v c is slightly reduced 
for the force of eq. [5]). This result is not surprising, 
since for r p ~ a the lattice potential effectively turns in- 
side out. A vortex segment now damps to a position that 
maximizes its distance from the nearest nuclei to the ex- 
tent possible against opposing tension forces. For r\ r = 
the vortex is trapped by the lattice as for the case of at- 
tractive nuclei. Since these results are independent of the 
value of rj r provided r\ r << 1, they are also independent 
of the assumption that r\ r is a constant. 

Recent calculations of E vn , using the local density ap- 
proximation |24] and mean field theory [25j are not in 
complete agreement. Though both calculations predict 
comparable interaction energies (~ 1 MeV) above an 
inner-crust density of ~ 10 13 g cm -3 , they do not agree 
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FIG. 3: Motion and pinning of the vortex in the x — z plane 
for the simulation of Fig. 2 with v s — 0. 
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on the densities where the interaction is the strongest. 
The inevitability of pinning is nevertheless a robust re- 
sult; it is insensitive to the exact density range in which 
E vn is strongest, the length scale of the potential, or even 
the sign of the potential. Within these uncertainties, the 
critical velocity above which vortices will translate is 



,1/2. 



p s na 



10" 



10 7 cm s _1 . 



(7) 



The lattice will trap the vortex for much smaller values 
of E vn as well, though v c will be lowered according to eq. 

Conclusions and implications. — Pinning of vortices to 
the lattice of the NS inner crust below a critical value 
of v s appears to be inevitable, whether the interaction 
is attractive or repulsive (and provided E vn exceeds the 
stellar temperature, typically ~ 0.01 MeV). For interac- 
tions of order ~ 1 MeV per nucleus, the pinned vortex 
lattice becomes unstable for v s > v c ~ 10 6 — 10 7 cm 
s , and vortices will translate approximately with the 
SF if the drag is low. For critical velocities this large, the 
pinned superfluid can store enough angular momentum 
to drive the giant glitches seen in pulsars (27] . If pinning 
in regions of the inner crust does occur, however, inter- 
pretation of the putative precession seen in some pulsars 
becomes problematic; as noted by Shaham [3], if even a 
small portion of the inner crust vorticity is tightly cou- 
pled to the solid, the star will precess much faster than 
the period of ~ 1 yr indicated by observations [7J. In 
the outer core, where the protons are predicted to form 
a type-II superconductor, vortices are expected to pin 
against a disorganized system of flux tubes and a sim- 
ilar difficulty arises in explaining long-period precession 
[5]. Anywhere there is pinning, however, vortices can 
creep (y << v s ) through thermal activation or quantum 
tunneling processes [3J |5S] not considered here. Vortex 
creep is incompatible with long-period precession if it is 
a high-drag process [T] [20] . 

The core neutron-charge mixture could be unstable 
to the formation of SF turbulence in a precessing star 
[T51 |2"5] , though magnetic stresses might suppress such 
an instability [3D]. [Turbulent instabilities might not oc- 
cur in the single-component SF of the inner crust] . These 
analyses used the drag description of vortex motion which 
does not correctly describe vortex motion through a po- 
tential; it is necessary to include non-dissipative pinning 
forces. The mutual friction force of eq. [TJ should be 
replaced with the total force per unit volume exerted on 
the SF by the vortex array. For a single-component SF, 
this force is (see, e.g., pQ): 



F = p s {v(v s ) - v s } x u, 



(8) 



where v(v B ) is local vortex velocity with respect to the 
background, averaged over many vortices, in response to 
a local flow v s . In the approximations made here, v = 
for v s < v c , and F reduces to the Magnus force per unit 



volume on pinned vortices. At finite v Sl however, vortices 
can move slowly through vortex creep, a process which 
dissipates energy at a rate p s v ■ (v s x oj) per unit volume; 
the dissipative force in the fluid would be small if v(v„) 
is nearly orthogonal to w s x u>, while the non-dissipative 
force would be, in any case, nearly equal to the Magnus 
force on perfectly pinned vortices. It would be interesting 
to study NS star modes under the assumption that creep 
proceeds with little dissipation, the opposite limit of that 
studied so far. 

I thank I. Wasserman and CM. Riedel for useful dis- 
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